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Abstract 

The inclusive production rate of mesons from the decay of b hadrons produced in Z° 
decays was measured to be Br(b— >0X) = 0.0282 ± 0.0013 (stat.) ± 0.0019 (syst.), using data 
collected by the OPAL detector at LEP. 
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1 Introduction 

In Z° decays, mesons are produced from three sources: direct production in the fragmentation 
process, decay of charm hadrons in TP^cc events and from the decay of bottom hadrons. 
According to the spectator model, Bu and Bj mesons decay to a meson via the cascade 
decay, b— >-c— Previous measurements of the branching ratio Br(Bu/ Bj— i>0X) were made by 
the CLEO |I| and ARGUS [0 collaborations and were averaged by the Particle Data Group to 
give the result Br(Bu/Bd^0X) = 0.035±0.007 [|. 

At LEP, where all b hadron species are produced, another mechanism exists where a Bg 
meson also decays directly into a meson. A measurement of the production rate in inclusive 
b hadrons, and in particular the rates of the individual b hadron species, can be used as an 
input to improve our knowledge of the Bg oscillation frequency. 

In this paper, we describe a measurement of the production rate of mesons from b hadrons. 
This includes both direct and cascade b hadron decays. mesons were identified in their charged 
kaon decay mode, and those not coming from Z°— i>bb events were suppressed by utilizing the 
relatively long lifetime of the b hadrons. 

2 Hadronic Event Selection and Simulation 

The data used for this analysis were collected by the OPAL detector from e"'"e~ coUisons at 
LEP between 1990 and 1995 running at centre-of-mass energies in the vicinity of the TP peak. 
Hadronic Z° decays were selected using the number of tracks and the visible energy in each 
event as described in reference [0. This selection yielded 4.41 million hadronic events. 

For optimisation of the selection of events and for some of the studies of systematic uncer- 
tainties, we generated 6.5 million 5-flavour hadronic TP decays (referred to as qq Monte Carlo). 
Monte Carlo events were also used to determine the selection efficiency for signal and back- 
ground. For determination of signal efficiency, we generated 80,000 Z°— >-bb events in which at 
least one of the b hadrons contained a meson in its decay chain. For background efficiency, 
we generated 1 million Z°— >^cc events. 

All these samples were generated with the JETSET 7.4 Monte Carlo program |^ with 
parameters tuned to the OPAL data 0. The heavy quark fragmentation was parameterised 
by the fragmentation function of Peterson et al. all samples were processed with the full 
OPAL detector simulation package and subjected to the same reconstruction and analysis 
algorithms as the data. 
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3 Analysis Procedure 



In each event, tracks and electromagnetic clusters that were not associated with a track were 
combined into jets, using the invariant mass algorithm with the EO recombination scheme and 
the parameter ?/cut set to 0.04 [fH]. The primary vertex of the event was reconstructed using 



the tracks and the knowledge of the position and spread of the e"'"e~ collision point [ pJ| . 

We searched for the decay 0^K+K~ in the hadronic event sample by combining two tracks 
to form a meson candidate. All two track combinations from the same jet with opposite 
charges were considered. Each of the two tracks was required to have a momentum of at least 2 
GeV/c. Due to the small opening angle of the two kaons, this momentum criterion is equivalent 
to a cut on the (j) momentum of about 4 GeV/c, which significantly reduces background from 
fragmentation while only accepting tracks in a momentum range with good separation between 
kaons and pions using ionisation energy loss (di?/dx) information [|T2|. The dE/dx for a kaon 
candidate track was required to be consistent with that expected for a kaon with a probability of 



more than 10% as defined in |T2|. To further reduce the pion background, the probability for a 
pion hypothesis was required to be less than 10%. These dE/dx selection criteria retain 51% of (j) 
mesons in the momentum range indicated above, while reducing the combinatorial background 
by 97%. Figure a) shows the momentum distribution of mesons from fragmentation, b 
hadron and c hadron decays in generator level simulation. The fragmentation meson spectrum 
peaks at low momentum but also has a tail to high momentum above the signal mesons. We 
therefore required the momentum of the mesons to be smaller than 25 GeV/ c. 

Since at this point, the hadronic data sample consisted mostly of non-bb events, we sup- 
pressed these events by means of a b-tagging algorithm, based on reconstructed displaced 
secondary vertices. An artificial neural network with inputs based on decay length significance, 
vertex multiplicity and invariant mass information |1^ , was used to select vertices with a high 
probability of coming from b hadron decays. Events were accepted if at least one of the jets 
was tagged by the neural network. The b-tagging selection was found to be 68% efficient for 
the signal, while rejecting 80% of the remaining background. 

The two kaon candidate tracks were fitted to a common vertex in three dimensions and the 
decay distance, the distance from the primary vertex to the reconstructed secondary vertex, was 
calculated. Figure |1| b) compares the decay distance distribution of mesons from fragmenta- 
tion and from b and c decays which passed the above selection criteria. In this plot, negative 
decay distance represents candidates with a reconstructed secondary vertex in the hemisphere 
opposite to the candidate's jet. Motivated by this plot, we rejected candidates with a negative 
decay distance. This selection left 62% of the remaining background events, while keeping 81% 
of the signal events. 

Candidates were accepted if their invariant mass was in the region: 1.011 GeV/c^ < Mkk 
< 1.027 GeV/ c^, which corresponds to about two standard deviations of the reconstructed mass 
distribution around the nominal mass. Figure |^ a) shows the invariant mass distribution of 
the K"'"K~ candidates. A total of 4297 candidates were found in the signal mass region. 



3.1 Background Estimation 

Two of the three sources of mesons mentioned in section are considered as background, in 
addition to the non- background where the candidate is formed from a random combination 
of tracks not both from a 0. To estimate the combinatorial background, we fitted the Mkk 
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Figure 1: Monte Carlo distributions of (a) true momentum and (b) reconstructed decay distance 
of mesons from b decays (solid), fragmentation (dashed) and charm decays (dotted). 



distribution to the form 

/(Mk) = (Mkk - 2 ■ Mk)" -{b + c-M^K + d- M^^), (1) 

where Mk is the mass of the charged kaon. The parameters a, b, c, and d were determined from 
a fit to the combinatorial background in the qq Monte Carlo sample, and then the function 
shape was normalised to the data in the mass region 1.060 GeV/c^ < Mkk < 1.225 GeV/c^. 
The background shape is superimposed on Figure |^ a). To verify that this technique estimates 
correctly the background shape, we compared the background-subtracted Monte Carlo sam- 
ple with the true distribution of Monte Carlo mesons (figure ^ b). We also compared the 
background-subtracted (p mesons in the data with the true Monte Carlo (p distribution (figure 
^ c). Both sets of histograms agree with each other very well, thus justifying the procedure for 
combinatorial background estimation. Subtracting the combinatorial background, we obtained 
2667 ± 40 mesons in the signal mass range, where the error arises from the statistical uncer- 
tainty on the combinatorial background estimation. The above background estimation ignores 
a possible small contribution from fo(980). As such a contribution is not well established, we 
consider it as a systematic uncertainty as described in section WA. 



The number above contains a small number of mesons from Z"— >cc decays and from 
fragmentation. The former was estimated using known and calculated branching ratios of 



Du, Dd and Dg into mesons, as in [jl4[. This calculation is based on isospin symmetry ratios 
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Figure 2: a) Invariant mass distribution of candidates passing all the selection criteria. The 
solid line represents the combinatorial background estimation as described in the text. The 
arrows show the signal and fit normalisation regions, b) Comparison of the distribution of true 
Monte Carlo (p mesons with that obtained by subtracting the combinatorial background shape 
from a fit to the combinatorial background, c) Comparison of the distribution of true Monte 
Carlo (f) mesons with that obtained in the data after background-subtraction. 
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between measured and theoretical branching fractions and was updated to include recent results 
§. We used: 



Br(Du^0X) = 0.0182 ±0.0025, (2) 
Br(Dd-^0X) = 0.0182 ±0.0025, (3) 
Br(D, ^0X) = 0.175 ±0.025. (4) 

Using the partial widths of Z° into cc and /d^ ^ ^, the branching fractions of charm quark into the 
different D mesons as measured by LEP |TB[ , we estimated the number of mesons which were 
produced from charm quarks to be 52746 ± 8322, where the error includes the uncertainty on 
all of the above numbers. With an overall efficiency to survive the selection criteria estimated 
at 1.48% and using Br(0^K+K-) = 0.491±0.008 [|, we obtained an estimate of 389 ± 61 
mesons from this source. The error, obtained from the above numbers, reflects the systematic 
uncertainty on this estimation. 

The fragmentation background was estimated directly from the qq Monte Carlo which 
had been tuned to reproduce the meson yield in the data [0]. It had also been shown |[TB|, 
that the shape of the Monte Carlo distribution agrees with that of the data, in particular in 
the region below 35 GeV. We estimate the number of mesons from this source to be 485±20, 
where the error is from the statistical uncertainty on the Monte Carlo sample used for this 
estimation. 



4 Result and Systematic Uncertainties 

The number of mesons which originated from the decay of a b hadron was calculated by 
subtracting the estimated background contribution from the number of mesons seen in the 
data. The result is = 1793 ± 50, where the error is from the statistical uncertainty on the 
background level. To calculate the branching fraction, we used: 

^'"^ ^ = 2.Wr,^.iV,.,...Br(^^K.K-)' 

where A^had is the number of hadronic events and e is the selection efficiency for mesons which 
originated from b hadron decay. With 4.41 million hadronic events, r^^^/Fhad = 0.21642 ± 



0.00073 PI and e = 0.068 ± 0.001, we obtained 

Br(b ^ 0X) = 0.0282 ± 0.0013 ± 0.0019, 

where the first error is statistical and the second is a 6.6% systematic error, obtained as discussed 
in detail below. 



4.1 Systematic Uncertainties 

Systematic uncertainties arise from the limited accuracy with which the branching ratios used in 
this analysis are known, from the uncertainty in the simulation used to determine the efficiency 
and from the background estimation. All uncertainties quoted below are relative, i.e. with 
respect to the measured value of Br(b— »0X). Table 1 summarises the systematic uncertainties 
detailed below. 
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Systematic Source 


5(Br(b -> 0X))/Br(b ^ 0X) (%) 


Branching fractions 


3.6 


dE/dx 


3.1 


Fragmentation cf) 


2.8 


Heavy quark fragmentation modelling 


2.4 


Combinatorial background 


1.6 


Monte Carlo statistics 


1.4 


b lifetime and decay multiplicity 


1.2 


Detector modelling 


1.2 


Total 


6.6 



Table 1: Systematic uncertainties 



Branching fractions 

The uncertainty on the known branching fractions of ^ K+K~, T^i^/T\^s.Ai Tcc/rhad, /ouda 
and Du d,s 'P resulted in a 3.6% uncertainty on the branching ratio of b^0X. As Bg mesons 
also decay directly to mesons, the efficiency for detecting mesons from the decay of a Bg 
meson is slightly different from that of Bu d- Since only a small fraction of Bg— »-0X is expected 
to come from direct decays, we increased the fraction of cascade decays in the Monte Carlo 
from the current value of 85% to 100%, and recalculated the efficiency. We also varied the 
production rate of b hadrons according to the known values [|15| . We obtained an additional 
uncertainty on the efficiency of 0.3%. 



Modelling of dE/dx 

To estimate the uncertainty arising from the modelling of the dE/dx selection criteria, we 
compared the efficiency of the dE/dx cuts in Monte Carlo events and in data. We reconstructed 
mesons without the dE/dx criteria in data and in Monte Carlo and compared the yield 
to that obtained in section |^. We obtained an efficiency associated with the dE/dx selection 
criteria of 49.5% in data and 51.0% in Monte Carlo. Hence, the relative uncertainty on the 
signal efficiency associated with the dE/dx cuts is estimated at 3.1%. A similar technique was 



used recently by OPAL |17|, and consistent uncertainties obtained. 



Fragmentation mesons 

The largest dependence on the Monte Carlo in estimating the background levels is related 
to fragmentation mesons. Our selection criteria were designed to reject fragmentation 
mesons, thereby limiting the possible uncertainty associated with their simulation. Since we 
are estimating the fragmentation background directly from the Monte Carlo, we would like 
to establish the level of accuracy of the simulation of this source. We took the small difference 
of meson yield between the Monte Carlo and the measured value [Q, and obtained a relative 
uncertainty of 2.7%. 

A possible shape difference will affect the selection efficiency due to the momentum cuts. 
We checked the effect of varying the Lund symmetric fragmentation parameters a (0.11-0.32 
GeV~^), and h (0.48-0.58 GeV~^). These changes produced a negligible effect. We also varied 
the parton shower A value in the range 0.13 to 0.31 GeV. This resulted in a relative uncertainty 
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Figure 3: Momentum distributions of (a) data (j) mesons after background-subtraction (dots) 
and bb Monte Carlo (histogram), and (b) data and Monte Carlo candidates from events that 
were rejected by the b-tagger or by the decay distance cut. 



on the fragmentation (p production rate of 1.6%. We also compared the shape obtained when 



using the HERWIG event generator [|T9[, and obtained a relative uncertainty of 2.4%. The 
uncertainty on the production rate of mesons from the decay of b hadrons associated with 
the estimation of the fragmentation background was estimated at 2.8%. 

To establish the agreement between the data and the Monte Carlo we show in figure ^ a) 
the momentum distribution of background-subtracted data and that of simulated mesons 
from the decay of b hadrons normalised to the data area. We also compared (figure |^ b) 
the momentum distribution of samples which came mostly from fragmentation by selecting 
events that were rejected by the b-tagger or by the decay distance cut. The data are shown 
after combinatorial background-subtraction and the Monte Carlo histograms are normalised to 
the same number of events passing the hadronic event selection. According to the Monte Carlo, 
80% of these mesons are from fragmentation. 

Heavy quark fragmentation modelling 

Heavy quark fragmentation was simulated using the function of Peterson et al. The heavy 
quark fragmentation model parameters were varied to change the mean scaled energy of weakly- 
decaying bottom and charm hadrons within the experimental range: {xE)h = 0.702 ±0.008 and 



{xe)c = 0.484 ± 0.008 [|T5|. These changes resulted in a 2.4% relative uncertainty. 
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Combinatorial background estimation uncertainty 

The uncertainty on the fitted shape parameters and on the normahsation gave an uncertainty on 
the background estimate. We also repeated the procedure with a function of the form /(Mk) = 
(Mkk — 2 ■ Mk)" ■ (6 + c ■ Mkk), (ci = in equation 1). This resulted in 35 additional events 
or 1.5% relative uncertainty. We also used the Monte Carlo shape directly (without fitting), 
by normalising the number of combinations in the mass range 1.06 < Mkk < 1-2 GeV/c^ 
to that found in the data. With this method, we obtained 1681 events (34 events more than 
our initial estimate). Since the above two estimates are similar, we took the larger of the 
two and combined it with the 0.5% of uncertainty from the fitting parameters, to obtain 1.6% 
uncertainty on the branching ratio associated with the background estimation. 

As the contribution from fo(980) was neglected in the background estimation, we assigned 
a conservatively large estimation of this contribution as a systematic uncertainty. We assumed 
that the width of the fo(980) is 100 MeV and that Br(fo(980)^KK) = 20%. We used an fo(980) 
multiplicity of 0.14 mesons per hadronic event and obtained the selection efficiency from 
the Monte Carlo. Combining these numbers, we estimate the contribution of fo(980) to the 
K+K~ invariant mass in the signal region to be 6 events. Assuming a 100% uncertainty in this 
estimate, this translates to 0.3% of relative uncertainty on the branching fraction of b^0X. 

b hadron lifetime and decay multiplicity 

The probability to reconstruct mesons from b hadron decays also depends on the efficiency 
to reconstruct secondary vertices in both hemispheres. This in turn is sensitive to the charged 
decay multiplicity and lifetime of the b hadrons. The Monte Carlo was reweighted to reflect the 
measured multiplicities and lifetimes |]T3[ . The uncertainty on these figures gave an uncertainty 
of 1.2% on the selection efficiency. 

Detector modelling 

The simulated resolutions of the tracking detectors were varied by ±10 % relative to the values 
that optimally describe the data, following the studies in The analysis was repeated and 
the efficiencies were estimated again. This source contributed an uncertainty of 1.2%. 

5 Summary and Discussion 

We have searched for direct and cascade decays of b hadrons into (j) mesons. The meson was 
reconstructed in its K"'"K~ decay mode and the fraction of Z°— i>bb events in the sample was 
enriched. We found the branching fraction for this process to be: 

Br(b 0X) = 0.0282 ± 0.0013 ± 0.0019. 

Of special interest is the production rate of mesons from Bg mesons. However, as the 
existing average of B"*^/" 0X measurements has 20% uncertainty, it is impossible to obtain 
a meaningful result for Br(Bs^0X). On the other hand, our result may be used to obtain an 
upper limit on the branching ratio of B^/° — 0X. We considered the extreme scenario in which 
Bg mesons decay to mesons only via Dg decays and used Br(Bs— >-DsX) = (90±33)% 0. We 
made a conservative assumption that b baryons do not contribute to the production rate and 
obtained an upper limit of: Br(B±/° 0X) < 0.0285 at the 90% confidence level. This limit 
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is consistent with the CLEO result (0.023±0.006±0.005) |l|], but is only marginally consistent 
with that of ARGUS (0.0390±0.0030±0.0035) g. 
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